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Under appropriate circumstances the electrons emitted from a superconducting tip can be en- 
tangled. We analyze these nonlocal correlations by studying the coincidences of the field-emitted 
electrons and show that electrons emitted in opposite directions violate Bell's inequality. We scruti- 
nize the interplay between the bosonic nature of Cooper pairs and the fermionic nature of electrons. 
We further discuss the feasibility of our analysis in the light of present experimental capabilities. 

PACS numbers: 03.65.Ud, 79.70.+q, 74.45. +c 
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Entanglement, at the heart of the foundations of quan- 
tum mechanics, has gained renewed attention with the 
birth of quantum information science [l[ . Here it is con- 
sidered to be a precious resource, as it is believed to be 
the key ingredient for the increased efficiency of quan- 
tum protocols as compared with their classical counter- 
parts. Finding sources of entangled particles is therefore 
of paramount importance. In quantum optics this is well 
known and routinely used. For example, in parametric 
down conversion 0] a pump laser beam incident on a non- 
linear crystal leads to the generation of entangled pho- 
ton pairs. In electronics, the field is much younger, but 
there are already a number of very interesting proposals 
to generate entangled electron states (see the reviews [1] ) . 
In this Rapid Communication, we propose an entangled 
electron source in vacuum, on the basis of a thorough 
analysis of the entanglement and the correlations of the 
electrons field-emitted from a superconductor. 

When a bias voltage is applied to a sharp piece of mate- 
rial, a strong electrostatic field is realized at the tip, caus- 
ing electron emission into vacuum. The ground state of 
superconductor is a fairly controllable macroscopic quan- 
tum state and provides a coherent and monochromatic 
electron beam via field emission from a superconduct- 
ing tip, as experimentally shown in Our analysis 
will show that electron coincidences in field emission can 
reveal electron non-local correlations due to pairing in 
the superconducting tip. Moreover, we shall see that, 
by orienting the detectors in opposite directions, one can 
optimize the fraction of entangled electrons in order to 
perform a test of Bell's inequality. 

Field emission thus enables one to study Bell's inequal- 
ity on electrons in vacuum by means of correlation mea- 
surements. Signatures of quantum statistics and of corre- 
lation can be unambiguously detected by coincident mea- 



surements. After the seminal result by Hanbury Brown 
and Twiss , bunching and antibunching of bosons and 
fermions have been measured in a series of important ex- 
periments @, 0, S, H, [m . Advanced technology in 
single-electron detection has made possible the observa- 
tion of antibunching in field emission [lo| . An experi- 
ment with a superconducting tip, similar to is what 
is needed to test our proposal. 

We will show that the spectrum of emitted electron- 
ics from a superconductor displays a remarkable inter- 
play between positive correlations (bunching), due to the 
bosonic nature of the Cooper pairs, and negative cor- 
relations (antibunching) due to the fermionic nature of 
the electrons that make up the Cooper pairs. An elec- 
tron pair exhibiting positive correlation is in a Bell state. 
The presence of a bunching-like behavior in quantum 
transport through multiterminal superconductor-normal 
metal structures has already been pointed out through 
an analysis of current noise [l^, [ij]. Superconductors 
were also shown to be a source of entangled electrons, 
whose generation and detection in hybrid conductors has 
been discussed in a number of articles [l^ . In particular, 
we draw attention to the work of Prada and Sols in [3] , 
where the angular distribution of the emitted electrons 
from a superconductor has been discussed. 

Our setup is sketched in Fig. [TJ The Hamiltonian in 
3D space reads H = Hs + Hy + Ht, with {h = 1) 
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FIG. 1: (Color online) Field emission of electrons and coinci- 
dent detections. The size of the emitting region at the tip is 
w and two of the emitted electrons are detected at ri and r2- 



where Hs is the Hamiltonian of the superconducting 
emitter, aks the fermionic operators of the quasiparticles, 
Hv describes the propagation of the electrons in vacuum, 
Cps are fermionic operators, and the (interaction) Hamil- 
tonian Ht describes the emission of the electrons from 
the superconductor into vacuum [l^i EBj [l3 ■ The energy 
of an electron in vacuum, e^, as well as that of a quasipar- 
ticle in the superconductor, w^, are both measured from 
the Fermi level of the superconductor, ^. A is the gap of 
the superconductor and the quasiparticlc operators a^s 
are related to the electron operators a^s by a Bogoliubov 
transformation [l8l | . The Coulomb interaction among the 
emitted electrons can be safely neglected since it becomes 
relevant at much larger current densities, as compared to 
those typical of the experiments relevant for the present 
work [l9|. 

For the tunneling matrix elements we take 



Tpk = h{p)g{p - k), 



gip) = (27r) 
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where g{p) characterizes the emitting region of size w 
and /i(p) the tunneling probability, that decays expo- 
nentially as the incident energy decreases [i^] , Ec being 
an energy scale that characterizes the energy dependence 
of the tunneling matrix elements. 

The emission process is dynamically described in 3D 
space with the Hamiltonian (IT|)-(I3]) and after a transient 
period it reaches a nonequilibrium steady state [2lj , with 
a stationary beam of electrons emitted from the super- 
conductor. Electron correlations are present in this beam 
17|. When the detectors do not resolve the spin states 
of the electrons, the probability of two joint detections, 
at (ri,ii) = (1) and {r2M) = (2), with ±2 > is pro- 
portional to 

/2)(2;1)= (V'1,(1)^1,(2)V'.,(2)V'.,(1)) 

Sl,S2 = T,i 

= 47(2; 2)7(1; l)-2|7(2;l)|2 + 2|x(2;l)p, 

(6) 

where ipsif,t) is the field operator of the electrons in 
vacuum, 

7(2;l) = (^|(2)^T(l)> = (^I(2)^i(l)>, (7) 
X(2; 1) = (^T(2)^i(l)> = -(^i(2)^t(1)>- («) 

The correlation function 7 describes the state of single 
electrons, and, in particular, the (spin-summed) one- 
particle distribution of the emitted electrons is given by 
p'^^-'(r,t) = 2^{r,t]r,t). The correlation function x de- 
scribes the emission of pairs of electrons (with opposite 
spins). A second-order calculation for the coincident de- 
tections at ti = t2 yields 



X(2;l) = / d^kukVk 
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where Uk and Vk arc the Bogoliubov amplitudes. Notice 
that, since UkVk = A/2uJk, X is proportional to the gap 
parameter A and vanishes when the emitter is in its nor- 
mal state. At zero temperature, there is no contribution 
from the quasiparticlc excitations and Eq. ([9]) is due to 
Andreev processes. In the absence of this contribution, 
the second term in ([6]) reduces the coincidence probabil- 
ity of finding two electrons close to each other within a 
small time delay, exhibiting antibunching. The pair cor- 
relation Xi on the other hand, enhances such coincidence 
probability. This_is_ relevant for the occurrence of posi- 
tive correlations 
in greater details. 
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Q{r,0) 



p(i)(2)p(i)(l) 
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at ti = t2, when the detectors are at the same distances 
ri = r2 = r from the tip, is plotted in Fig. [2] as a func- 
tion of the angle 6 between ri and r2, for normal and 
superconducting emitters. Here, kp = ^J2ra[i = 2'k/\p. 
The effects of superconductivity are manifest: a bunch- 
ing peak appears at ~ tt. Its origin is clear from the 
expression of the Andreev process ([9|) (7 is negligibly 
small at 9 ~ tt). This shows that electrons with opposite 
momenta k and —k are emitted in pair through a vir- 
tual process and propagate with momenta pi and p2 in 
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vacuum, in approximately opposite directions (with un- 
avoidable diffraction effects governed by the size of the 
emitting region w). The couple k and — fc reflects the 
Cooper-pair correlation in the emitter. Notice that the 
integrand of x in ([5]) is symmetric under the exchange 
k ^ —k. This is because the Cooper pair is in a sin- 
glet spin state. This symmetry yields bunching, which is 
observed in opposite directions. 

Bunching is therefore a signature of excess singlet 
pairs, when the emissions take place in opposite direc- 
tions. It is then of great interest to discuss the nonlocal 
aspects of the phenomenon 0, ■ The spin state g of 
the pair of emitted electrons is 



(7(2; 2)7(1; 1)-|7(2;1)P) 1 
+ 2(|7(2;1)P + |X(2;1)P) 



(11) 



where |*-) = (|ti) - |iT>)/V2 is the singlet state [the 
normalization factor is given by the two-particle distri- 
bution p(^'(2;l) in ©]. Therefore, in general the en- 
tanglement of the singlet component is masked by the 
background. 

The degree of entanglement is related to the height 
SQ = Q{r, 6*) — 1 of the bunching peak at = tt, which, 
for kpr » 1, » Ec, and r/kpw"^ ^ 1, is given by 



32Kf{\A\/Ec) 



H, 



(2) 



IW 
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where K,j{z) is the modified Bessel function of the sec- 
(2) 

ond kind. Hi (z) is the Hankel function of the second 
kind, ^ = fcj?/7rm|A| is Pippard's length, characterizing 
the correlation length of the superconductor, and A is 
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FIG. 2: (Color online) Normalized coincidences Q{r,9) vs 6, 
for normal A = and superconducting A 7^ emitters. Ob- 
serve the large bunching peak at ~ tt. 



a smooth bounded function of w, such that A ~ 1 for 
ui > Af- The higher the bunching peak, the larger the 
entanglement, and eventually Bell's inequality can be vi- 
olated. Notice that the electrons of each pair are emitted 
in opposite directions and one need not argue how to 
separate them. 

Figure [3] displays the behavior of the bunching peak, 
by scrutinizing the role played by the parameters describ- 
ing the system. Clearly, the value of A is very significant 
for entanglement, as the effects of superconductivity are 
enhanced. By increasing |A|, the gap becomes wider, 
and entanglement is enhanced. The parameter Ec ap- 
pearing in (O works like a filter: by decreasing Ec, the 
contribution of single-particle emission is suppressed, the 
background is reduced, pair emissions become dominant 
and entanglement is enhanced. 

The effects of w on entanglement are interesting to 
discuss. Electron pairs emitted from a smaller region 
bunch better and are more entangled. If the emitting 
region is larger, there is less guarantee that coincidence 
electrons originate from a common Cooper pair, and as a 
consequence entanglement is reduced. This explains the 
role of the ratio between the size of the emitting region 
and the extension of a Cooper pair w/^, appearing in 
the formula for the bunching peak p2p and governing 
the entanglement of the emitted pairs. 

Finally, let us focus on the effects of propagation. A 
smaller value of r yields more entanglement. This is be- 
cause the wave packets of the emitted electrons spread 
as they propagate. Even if two electrons are detected at 
the same distance in opposite directions, this does not 
ensure that the two electrons originate from a common 
Cooper pair: there is an ambiguity to the extent of the 
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FIG. 3: (Color online) Peak value of Q(r, 6) a.t 6 = tt. The 
line at Q = 3/2 indicates the entanglement threshold (above 
which the pair of electrons is entangled) and the one at Q = 
\/2/ (\/2 — 1) — 3.41 is that for the violation of Bell's inequality 
(above which Bell's inequality is violated). The parameters 
are the same as in Fig. O and ^ ~ 33.8 \f- 
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spreads of the wave packets. Due to free-space propaga- 
tion, the uncertainty at time t ~ mr/kp is < \ftjrn and 
this value should be smaller than ^ for the two electrons 
to bunch. The bunching peak (fT2|) actually decays like 
~ kp^ jr for r kp^'^ (with oscillation), but the length 
scale fc_F^^ is much longer than the extension of a Cooper 
pair ^, and the slow decay reflects the divergence in 
the quasiparticle spectrum. The oscillations of 6Q (be- 
low the entanglement threshold) shown in the last panel 
of Fig. [3] for large values of r, are due to the asymptotic 
behavior of the Hankel function Hq ' . 

It is important to check to which extent our results are 
robust in a non-ideal situation. To this end we analyzed 
both static fluctuations of the diameter w and the posi- 
tion ro of the emitting tip. In particular, fluctuations are 
important only when they become of order of ^. More- 
over, one can show that the angular dependence of the 
peak is oc exp{— Sfcl^w^ sin^[(7r — 0)/4]}. Therefore, the 
effect should be visible as far as kpw56 < 1, where S9 
is the angle deviation from tt in the emission of the two 
correlated electrons due to local imperfections of the tip. 
This implies a maximum tolerable value of the roughness 
of the order of 1/Sk = l/kpS9 ~ w. 

In conclusion we have shown that field emission from 
a superconducting tip can provide a source of entangled 
electrons in vacuum. Besides being of great importance 
for the generation of entanglement in electronics, we be- 
lieve that a test of Bell's inequality on field-emitted elec- 
trons is of interest in itself. Moreover, this would be a 
remarkable example in which the interplay between the 
bosonic nature of Cooper pairs and the fermionic na- 
ture of electrons is brought to light. Although all the 
ingredients to experimentally observe our predictions are 
already available, our analysis shows that stringent re- 
quirements should be met, as suggested by Fig.[3l A large 
energy resolution and a tip material with a large value 
of the gap are certainly desirable. Also, energy selection 
close to the Fermi level would enhance correlations. 
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